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Abstract—The analytical model previously developed by 
the authors for backscattering from a random air- 
dielectric interface is used to investigate the effects of 
coherence in the backscattered field on the monostatic 
radar cross-section (RCS) of the sea. The sea-surface is 
represented by several empirical wave-number spectra 
for a fully developed sea, namely, the Pierson-Neumann 
Spectrum, Pierson-Moskowitz Spectrum, and_ the 
Elfouhaily Unified Spectrum. A Generalized Normalized 
Radar Cross-Section (GNRCS) is introduced as a metric 
of the combined effects of coherence and antenna 
radiation characteristics. By employing a small-slope 
approximation (SSA) of the sea surface, the relative 
powers of the signal returns are computed and compared 
to the signal-power returns based on data-derived 
normalized RCS X-band (NRCS) (o°) values for the 
corresponding simulated sea state. The objective of this 
analysis is to assess the effects of spatial coherence on 
backscattering using a GNRCS through validation with 
o° obtained from measurements of sea backscatter. 


Index Terms— Generalized Normalized Radar Cross 
Section (GNRCS), Antenna Reciprocity, Ocean 
Backscattering, Antenna Radiation Pattern. 


I. INTRODUCTION 


Normally, the assessment of radar-clutter returns from a 
distributed target (that is land, sea, etc.) is based on models 
that discretize the scattering contributions from within the 
antenna footprint [1]. The surface within the main beam 
footprint is treated as a random distribution of point 
scatterers with no mutual interaction. The individual 
contributions from the point scatterers are taken as 
proportional to the respective incident signal powers and 
summed. To maintain consistency, the vertical displacement 
of the surface is modeled as a random process, with no 
consideration for the spatial correlation within the beam 
footprint. Based on these assumptions, the relative power of 
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the radar clutter is estimated by evaluating the surface 
integral of the gain-squared divided by the fourth power of 
the range to the local region encompassed by the antenna 
footprint, and multiplying the result by the normalized radar 
cross-section (NRCS) [1], or o°as it is commonly known. 
Thus, stated explicitly, 


2 G7 
Pe Z A 5 f (8,4) ods. 
Pans (47) R 


The NCRS is the most important metric in the 
characterization of radar clutter and its properties are 
assumed to be independent of the antenna radiated-field 
features within the illuminated region of the surface. For the 
most part, o° is considered to be separable from the variables 
entering into the surface integral and is determined by its 
direct inversion. 


Obtaining o° with this technique requires that the 
illumination within the area contributing to backscattering be 
constant. However, when the clutter originates from a non- 
uniformly illuminated patch of the sea surface (say the entire 
antenna footprint) the validity of this inversion technique is 
questionable since it entirely ignores any spatial correlation 
among the clutter returns from different regions of the 
antenna footprint. In addition, the assumption associated 
with the inversion technique fails to take into account of the 
statistical and physical properties of the surface as well as of 
the antenna radiation pattern. This is represented in a number 
of examples in the literature. For instance, in [2], [3], [4], and 
[5] the gain and polarization were assumed to be uniform 
throughout the antenna footprint. Here the gain and 
polarization contributions from the center of the main lobe of 
the antenna pattern are taken to be the only values of interest. 
In [6], [7], and as discussed in [8], the uncorrelated scattered 
field assumption is enforced based on the decorrelation times 
of the radar returns. 


It has been demonstrated that the measured radar cross- 


section (RCS) per unit area of a rough surface depends on 
both the kinematic and stochastic properties of the rough 
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surface as well as on the properties of the receiving 
(transmitting) antenna radiation and polarization patterns [9]. 
In addition, formulations of o°, which incorporate the 
effects of spatial coherence, have been investigated in [10] 
and [11]. A most important feature is illustrated in [10], 
where the author distinguishes the difference between the 
NRCS generated by a plane wave and an antenna 
illuminating a rough surface. Herein, the author establishes a 
criterion for the NRCS being equal for both types of 
illumination; provided that the radii of curvature of surface 
are small compared to the antenna footprint. However, these 
attempts assume an incoherent scattering model and/or only 
consider the radar returns from the beam center of the 
antenna pattern. An approach to modeling backscattering 
from randomly perturbed surfaces is presented in [8] and is 
capable of incorporating these effects. That model was 
subsequently applied in [12] in a study of spatial-coherence 
effects of the electromagnetic field backscattered by an 
irregular dielectric-air interface. In this paper, an extension 
of the formulation of [12] is applied to the analysis of radar 
ocean clutter using a modified version of the small-slope 
approximation (SSA). The wave-number spectra of the ocean 
surface are represented by the Pierson-Moskowitz, the 
Pierson-Neumann, and the Elfouhaily unified wave-number 
spectra for a fully developed sea. As a measure of back- 
scattered power, o° is replaced by the Generalized NCRS 
(GNRCS) per unit area which incorporates the effects of 
correlation of the backscattered field as well as its 
dependence on the relevant antenna parameters as observed 
within the footprint. Stated explicitly, embedded in the 
formulation of the GNRCS are: 


e the statistical and physical properties of the surface; 

e the kinematic description of the surface; 

e the gain and polarization generated by the antenna. 

Numerical results are presented in terms of the average 
relative received power as a function of wind speed, wind 


direction, antenna radiation pattern, and polarization, and 
compared with the corresponding relative received powers 
based on measured o° at X-Band. Lastly, we establish the 
product of the antenna gain and GNRCS as a metric for the 
combined effects of spatial coherence of the surface and 
antenna radiation characteristics within the footprint. The 
objective is to assess the effects of spatial coherence on 
backscattering using the GNRCS through validation with 


0 . 
o obtained from measurements of sea backscatter. 


II. RCS BASED ON AN INCOHERENT SCATTERING MODEL 


An antenna at mean height H above the ocean surface is 
shown in Fig. 1. The Cartesian coordinate system serving as 
a reference for the spherical coordinates 6 and ø is 


indicated by the red axes. The antenna can be beam steered 
in elevation to @=6, and in azimuth ø = ø .The antenna 
gain will have its maximum atO=0, ¢= $ corresponding 


to the position vector R . The position vector R is directed 
from the phase center of the antenna to a given point on the 


rough surface with magnitude R = 4/x° + y° +(H -ç(x, y)? . 
Forç(x,y)=0, R=R, =4x +y +H’. 


Circular Dish Antenna Uniformly 
Illuminated on Transmission 


Rough Surface Modeled as a Random 
Process of the Horizontal Coordinates 


H - distance between Antenna and 
Mean Level of Surface 


Fig. 1. Scattering geometry and antenna footprint 


The solid angle Q determines the area of the beam 
footprint A on the surface. A beam-unresolved scatterer is 


defined as a scatterer that subtends a solid angle Q 


Antenna 


scatterer 


much smaller than Q . As a result, the field incident on 


Antenna 
the scatterer can be approximated by a plane wave and the 
receiver-to-transmitter transfer function can be expressed by 


[8]: 
be _ 20 j G(9,9) oO -2 jkyR 
x 2 P, ‘P,e , 
aans ky 4aR 4r (1) 
where 4,,,,; iS the stimulus from the transmitter, b,,. is the 


response due to the impinging plane wave, the scattering 
cross section of a scatterer o is approximated by the small 


patch of the illuminated surface and G(6,¢) is the available 
antenna gain. The p, and p, are the (complex) polarization 
vectors, respectively, of the transmitted and scattered fields, 
all evaluated at (0,¢). When Q <Q we speak 


of a beam-resolved scatterer for which the calculation of the 
receiver response requires an integration of (1) over Q 


Antenna scatterer ? 


Antenna * 
Radar clutter estimates are generally based on an incoherent 
integration (summation). This amounts to summing the 
squared magnitude of (1) weighted with the differential cross 


*) =o°dS, where dS is the differential 


section d (ofp. "P, 


area. As a result one obtains 


2 


(2) 
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which is just the average of the ratio of received-to- 
transmitter available power. 


The usual approach in finding o° from (2) is to assume 
that within the main beam (for all practical purposes at the 
beam peak) o° is constant and simply invert (2). One 
problem with this procedure is that it fails to take account of 
contributions originating outside the main beam however 
defined. Within the main beam, no account is taken of 
incoherent addition as a potential source of errors. It appears 
that the justification for incoherent processing rests mainly 
on its relative simplicity. Another common idealization is to 
entirely ignore the intrinsic dependence of o° on the antenna 
radiation field. This dependence has been identified in [8], 
[10], and [12], and is incorporated into the GNRCS 
discussed in the sequel. 


Ill. 


As given in [8] and elaborated on in [12], the receive- 
transmit transfer function linking the transmitter output 
aans With the response to the backscattered field b,,. can be 


trans 


COHERENT SCATTERING AND THE GNRCS 


put into the form: 


be _ -2 G(0, 
rec T japans Í ( g) p x S jise (9, $) P p,e 


—P, ~22koR g2R 
4zR 


(3) 


0“ trans A 


where Sas (0,9) isa 2x2 scattering matrix determined by 
the back-scattering characteristics of a rough air/dielectric 
interface. The scattering parameters (in the network sense) 


è 2 
Arans aNd b,,. are normalized such that |b, Ta 


represent, H the received power and the 
transmitter available power. The transmitter polarization 
vector p, is defined by 


2: 
and ja 


(4) 


2 
and |F a (8,9) represents radiated power per steradian. 


Accordingly, the available antenna gain is 


1 Fad (9, $) 


a 


trans 


G(6,¢)=4 


Note that the expression in (3) is the antenna reciprocity 
relationship for a given antenna irradiated by an arbitrary 
electromagnetic field. Equation (3) was derived from the 
now Classical treatment which utilizes network parameters to 
give a full description of the electromagnetic phenomena of 
any antenna. This includes parameters like gain, polarization, 
aspect ratio, radiation efficiency, scattering, etc. Further 
details regarding the derivation of (3) are addressed in [13], 
[14], and [15]. 


Based on (3), the average of the ratio of received-to- 
transmitter available power is 


2 


SOO), -Sas (0,9) p,e RPR 


(5) 


To assess the consequences of approximating coherent 
scattering by incoherent scattering, we transform (5) into a 
form resembling (2). To this end, we first represent the 
differential d*R in terms of its projection on the xy plane 


d°R = s(0,) dxdy = s(0,9)dS = |1+|V,¢(x,y)| as 


and define a new scattering variable 


Q(0,) = 5(9,0)B, Sas (8,0) p, e CA. 


Factoring (5) and averaging only the product of the 
scattering variables yields 


i 


Arrans 


)- aji ie all D(a (9,¢)Q(0',¢'))ds‘ds. 


After rearranging the terms in the first integrand, we get 


2 7 IJ Th ai jee 


ee ji ae o(S)ds. 


The last equation has the form of (2) with o° replaced by 


-EEA (9(6,6)0(619)) 8 
a (6) 


which we define as the GNRCS. The GNRCS is not only a 
function of the radiation field (including its polarization) of 
the antenna, but also of the kinematics and statistics of the 
surface. As can be seen from (6), in regions where G(0,¢) 


is constant it can be removed from the integrand and 
cancelled with the gain function in the denominator. In that 
case, the GNRCS becomes independent of the antenna gain. 
However, this cancellation applies only to the gain and not to 
the polarization. Thus, in general, the GNRCS is polarization 
dependent. In view of (5), this dependence enters as an 


Q(9,9)Q(4',9')) as" 
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interaction between the antenna polarization and the 
scattering characteristics of the surface. 


Note that (6) reduces to o° (0,¢) when 


(2(0,9)0(0,6"))~ 7]? (0.8) 5(o-2), 


which corresponds to incoherent scattering. That is, if all of 
the surface properties are known to be statistically 
uncorrelated, then 


o(S)=0(8,6) > 0° (0,6) = 427| (0.9) 


2 


(7) 


The parameter y is some measure of the root-mean 


squared (RMS) height, slope, and/or horizontal correlation 
distance of the surface compared to k, . On physical grounds 


f (0,9) must encompass the polarization of the equivalent 


fields and constitutive properties of the surface in A. From 
an empirical standpoint, the o° obtained from the ocean (or 
from any rough surface) is just the average or mean value of 
a(S) given by 


J Jo(s)as 


;A=| fas 
A, (8) 


In examining (6) or (7), it is seen that properties of the 
surfaces are generally not separable for direct inversion for 


computing o°. It is the opinion of the authors that the best 


way to characterize ocean clutter is through the received 
power at the port, because that is what is actually measured, 
and not o’. 


IV. SSA OF THE SCATTERING MATRIX AND THE GNRCS 


The preceding rather general approach to surface 
scattering follows closely that presented in [8] (see Appendix 
A). These results are purely formal for they omit the crux of 
the problem, which is determining the scattering matrix (see 
Appendix B). Again, as outlined in [8] in its full generality, 
we are dealing with a complex boundary-value problem so 
that analytical results invariably employ some form of 
approximation. This paper employs the simplification 
introduced in [12] using an approximation closely linked to 
what has been referred to as the Rayleigh hypothesis [16]. In 
our interpretation, it amounts to approximating the fields 
scattered by an irregular boundary by the Fourier transform 
of the field scattered by an equivalent planar surface. This 
approximation is difficult to justify theoretically so that its 
validation must come indirectly from experimental data. 
Using this formulation, one can find the scattering matrix 
using perturbations either for small amplitudes (that is, 


[kos <1)) or for small slopes, 


len Set 


V(x, y)|< tveles)= K ay 


Only the SSA, or rather a modified version of it, is used in 
this paper and follows closely the exposition in [12]. The 
conventional formulation, using the SSA for scattering by an 
irregular boundary, is presented in [16]. Unlike in the present 
paper, the emphasis is on scattering per se with the antenna 
playing only an indirect role. Under the modified SSA the 
scattering matrix, restated directly from [12], Eq. 24, is 


Sas (9,9) = 


k, cos 8 


Eo [ketoos] 0) 


where I'(@) and T”(0) are the Fresnel reflection 


coefficients for parallel and perpendicular polarizations, 
respectively, at the mean level of the surface H in Fig. 1. 


Note that unless lkas (x,y) <1 


, the scattering matrix is a 


nonlinear function of k,¢(p). This nonlinear dependence 
corresponds to large vertical surface displacements (wave 
heights) compared to the radar wavelength. It is worth noting 
that ç( p) is a relative displacement with respect to the mean 


height H and is not affected by a displacement of the 
irregular surface in the vertical direction in this treatment. 


Ve(xy)] in (6) to 
unity, and since the antenna is sufficiently far away from the 
surface R ~ R,—¢(p)cos@, (6) assumes the form 


Small slopes force s(0,¢) =,/1+ 


2 -j2kR, 
_ 4k, Roe 


a(S)= 


1G(0,¢) 
(eee (p, Saisie (8,9): p,e), ~ 
Re \ (Pr Si (0, 8)-pre eh”) F 


(10) 


The terms in the expectation operator in (10) are direct 
complex-conjugate pairs in which the products can be 
simplified to 


Pr Sase (0): P, = sin| kag (p)cos8 | f,(0,¢). (11) 


where f,(6,¢) is the polarization factor under the modified 


SSA, which is equal to 
f,(8,¢) =cos6| p; (8,9) (8)- p; (0,¢)0"(8) |. a 


The product of sin| kaç (p)cos8 | and es)? in the 


ensemble expectation operator in (10) is also represented by 
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a set of complex conjugate pairs. This product is expressed 


4 . 4 y 
as l (EAEEREN eret, 
Y, (0,9) =sin| kig (p) cos 0 |e. {k k k ka} 


Kiisi = eltos(p)cosd | 


= {3k, cos 8, 3k, cos Ø, k, cos 8, k, cos 8}; 


~-~ m~~ 


= {3k, cos 0',k, cos 6’, k, cos 0’, 3k, cos 0’. 
2j 03) (16) 


A simplification is achieved when the narrow-beam 
approximation is applied ec) with the result For the purposes of this investigation, the PDF of ¢( p) is 
R, =R+ p- Asin (14) chosen to be a Gaussian random process with autocorrelation 


function R(p- p) and y (x) wave-number spectrum. 


The unit vector A, is defined as @ = X, cosġ +y sinó , Equation (16) then assumes the form 
where é is the antenn sami beam Stents angle. If the (v (0,0) £: (0, ’)) ak ` a n 
appropriate substitutions are made in (11) using (12), (13), A mi (17) 
and (14), then the GNRCS under the SSA is expressed as 

ako f, (8,4) E 
S)= EROS: eed Fy, T i 
a(S) oe i a o  R(p-p)|| Ke, 
sin $ -=k R( = ) lop -k 
JJO (0g) PeR (w (0,9): (0',9")) as. ADETE A 
(is) = kyon. (e cos? 0 + B? cos? 0'— 2a, 8, cos 0 cos g Zee) 
Equation (15) relates the material properties of the surface (18) 


and the electromagnetic properties of the backscattered field 
within A, to the spatial coherence of the surface as 


where /@,,@,,@,,@,' =}3,3,1,1' and 
characterized by the quantities (E, (0,6) ¥(8.¢')) and { 12 2 Q3 x} { \ 


{Bis Bo» Bs» Bs} ={31,1,3}, and o% is the square of the 
. The term (Y, (0,0): (0',¢')) is the RMS vertical displacement ¢(p) . Substituting (17) and (18) 


expected value of the Hermitian product of the phase-factor into (15) yields 
Y ,(0,¢), which is related to the impedance seen at the 


art sinĝ( p-p )A 


antenna port due to the vertical displacement ¢(p). The (s) _ 4ko f. (0.9) aol jel (0.0) f(A. #)e ~ j2ky sind (p-p') 

name "phase-factor" was selected because ¢(p) enters into 1G(0,9) 

the argument of the complex exponential in (13). An 1 “ -cn ({ * cos’ 0+ B? cos’ 6 
expression for a (0,9) Y; (0',¢')) can, in general, not be PAS mee ~2ar, B, cos Ocos of Bee) dS 
given explicitly until the joint probability density function E (19) 


(PDF) of ¢(p) is specified. For the ocean surface, the PDF 


can be constructed empirically through tedious 
experimentation. The results of these methods have not yet of a Gaussian random surface illuminated by a narrow-beam 
produced reliable PDFs and are not employed in this antenna. Note that the integrand includes the factor 
investigation. go sin 8(p-p )A 


where Goms = Koo... Equation (19) represents the GNRCS 


, which, as shown in [12], in the limit of 
Continuing onward, (E, (0,6) ¥:(8.9')) is expressed constant gain (and polarization) over the antenna beam 


as the sum of the characteristic functions for the different footprint, leads to Bragg scattering. 


spatial spectral weights k, and k, , which are directly Recall that Eq. (19) is valid only for surfaces that obey the 


derived from the inner and outer algebraic product terms of Constraint 


V, (x,y) «1 at every point on the surface. 
Y,(0,¢) 2 (0,9) or, Thus, (19) does not apply to spatially uncorrelated surfaces 
since it diverges when R(p-p ) is replaced by (p-p). 
To further illustrate this, it is known that for uncorrelated 


surfaces the mean-square slope VRP), cannot be small. 


5 
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Hence, a(S) cannot be based on an incoherent scattering 
model. On the other hand, when in addition to small slopes 
one also assumes small amplitudes, that is with kos ( p) «1 


, the exponential in (19) can be replaced by a linear function 
of R(p-p) so that (19) becomes compatible with an 


incoherent scattering model. As shown in equation (36) of 
[12], this compatibility is a general feature of the small 
amplitude approximation (SAA), independent of the 
constraint on the surface slopes. This feature does not imply 
that the SAA requires incoherence. The choice of the 
correlation function is arbitrary and can be chosen to model 
surface scattering where coherence is essential, as, for 
example in Bragg scattering [12]. 


V. SIMULATION OF THE GNRCS AND COMPARISONS WITH 
PUBLISHED DATA 
b /a 


reč. trans 


Plots of ( 


a as a function of the (elevation) 


beam-steering angle computed using (5) and (6), with o (S) 


represented by (19), are shown in Figs. 2 and 3. The plots 
are for the three aforementioned wave-number spectra and 
are identified in the legend with an "S" for simulated. Also, 


b/a a (or (2) 


rec trans 
derived from published data of o° for 4 m/s wind speed at 
crosswind. The plots for the horizontal polarization are 
shown in Fig. 2 and those for the vertical polarization are 
shown in Fig. 3. Similar results for wind speeds of 5 m/s are 
shown in Figs. 4 and 5. These simulations were conducted 
for the given sea-surface conditions at X-Band (10 GHz) and 
dielectric constant £, = 55.85- j37.71. The antenna is a 


shown for comparison are plots of ( 


circular dish, uniformly illuminated on transmission, and 
with a 16-wavelength diameter (0.48 m) placed 1,000 
wavelengths (30 m) above the sea surface (Fresnel number 
0.064). The average relative received power at the antenna 
port was computed for a given beam-steering angle ranging 
from 0° to 75°in 1° increments in elevation with the 
azimuth steering angle fixed at 90°. The simulated curves 
for o(S) in each figure are plotted in solid lines, with 


magenta for the ocean surface modeled by y(x) for the 


Elfouhaily Unified spectrum [17], light green for the 
Pierson-Neumann spectrum [18], and blue for the Pierson 
Moskowitz spectrum [19]. Details of the corresponding 
wave-number spectra can be found in [18], [19], and [17]. 
The data representing the median RCS per unit area of the 
sea surface, o°, used for the comparison (referred to as 
baselines curves) were obtained from various published 
accounts [20], [21], and [22]. These median values of the 
NRCS were measured at wind speeds of 4 to 6 m/s for a 
given elevation antenna beam-steering angle for crosswind 
wind direction and are identified in the legend with an "E" 
for experimental in each figure below. 


b /a 


rec trans 


The data used to derive the ( 


; ) using Long, 


[20], are depicted in Figs. 2 and 3 with a solid-black line 
with squares. The curves generated from data published by 
Daley, et. al., [21], are shown with a dashed black line with 
circles and the curves generated from Plant, et. al, [22], are 
shown with a dashed dark-green line with inverted triangles. 
A similar format is employed in Figs. 4 and 5 with the 
addition of a different data set as provided by Long and 
depicted with a dashed black line with squares, and the data 
set generated from Plant for 6 m/s uses green lines with 
triangles. The data from Plant were used in Figs. 4 and 5 to 


b /a 


rec trans 


confirm the trends of the simulated ( 


2 
) curves. 


b/a 


rec trans 


The code used to compute ( 


2 
) for each wave- 


number spectrum was initially developed in MATLAB 
where a method of integration derived from the 2D circular 
convolution was used to determine the GNRCS under the 
SSA. The auto-correlation function of the surface 


R(p—p) for Figs. 2 through 5 was computed from a 2D 
fast-Fourier transform (FFT) of the ocean wave-number 
spectrum y (x) fora grid size of 2048x2048. To accelerate 


the calculations, the code was ported into FORTRAN 90, 
and the integration portion was coded using the Message 
Passing Interface (MPI) to take advantage of the highly 
parallelizable nature of the algorithm. The calculations were 
performed using Cray XE6's made available through the 
DoD Supercomputing Research Centers at the US Air Force 
Research Laboratory at Wright-Patterson AFB and the US 
Army Engineering Research and Development Center at 
Vicksburg, Mississippi. 

The simulation curves in Fig. 2 proved to be promising for 
horizontal polarization, showing as little as a 4-dB deviation 
in the average relative received power as compared with the 
associated baseline values for smaller elevation angles (0 ° to 
40°). For the larger elevation angles (50° to 75°), the 
deviations from the baseline relative power curve generated 
from Long, [20], for the Pierson-Neumann and Pierson- 
Moskowitz models were as large as 10 dB. The Elfouhaily 
model appears to provide the best overall performance with a 
maximum deviation of 5 dB established between baseline 
curves for Long, [20], and Daley, [2 1]. It should be noted 
that the baseline curve for Plant, [22], deviates from the 
remaining curves in Fig. 5. In this case, Plant's baseline 


b/a F 


rec trans 
values whereas the baseline values associated with the data 
from Long would represent an upper bound. 


would represent a lower bound on the o° or ( 
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Or ; 
= Elfouhaily Spectrum 
sS —— Pierson-Moskow itz Spectrum 
=~ Pierson-Neumann Spectrum 
-20/ | =- Plant Show ex X-Band Fixed U = 4 m/s hh (2005) | 
E —& 4FR--Daley 8910 MHz 4 m/s hh (1971) 
== Long 8910 GHz 4 m/s hh (2000) 
-40- i 


2, 
IPsec rane! > (4B) 


0 10 20 30 40 50 60 70 
Elevation Beam Steering Angle (Degrees) 


Fig. 2. Relative Received Power vs. Elevation Beam Steering Angle of a 
horizontally polarized circular dish above a fully developed sea surface with 
wind speed of 4 m/s 


Or ; 
= Elfouhaily Spectrum 
sS —— Pierson-Moskow itz Spectrum 
20- |— Pierson-Neumann Spectrum i 
" = Plant Show ex X-Band Fixed U = 4 m/s vv (2005) 
E =¢ 4FR--Daley 8910 MHz 4 m/s vv (1971) 
== Long 8910 GHz 4 m/s vv (2000) 
-40 - 1 


2. 
Prectrans! > (dB) 


o 10 20 30 40 50 60 70 
Elevation Beam Steering Angle (Degrees) 


Fig. 3. Relative Received Power vs. Elevation Beam Steering Angle of a 
vertically polarized circular dish above a fully developed sea surface with a 
wind speed of 4 m/s 


The results for vertical polarization shown in Fig. 3 are 
similar to those of Fig. 2, where the simulated 


(|b /a 


rec trans 
baseline curve from Long within the range of elevation angle 
of 0° to 40°. In re-examining Fig. 3, it is seen that there is a 
larger deviation from the baseline curve of Long as well, as 
much as 12 dB for all of the simulation curves. Again, the 


(|b /a 


rec trans 
appears to be the overall better estimate for the higher 
elevation angles, with a smallest 2-dB deviation compared to 
the lower bound and an 8-dB deviation compared to the 
upper bound. 


‘| curves also exhibit good agreement with the 


*) generated from the Elfouhaily spectrum 


The simulation curves in Fig. 4 proved to be acceptable for 
horizontal polarization, showing as much as a 7-dB deviation 
in the average relative received-power responses when 
compared with the associated baseline values for elevation 
angles of 0° to 40°. For elevation angles of 50° to 75°, the 


deviations were as small as 4 dB from the baseline curve 
generated from Long for the Pierson-Neumann and Pierson- 
Moskowitz models. The Elfouhaily Spectrum again provides 
the best overall performance with a minimum deviation of 2 
dB from the baseline curve for Long and Daley. Note, the 
baseline curve at 4 m/s for Plant is down significantly from 
the remaining curves in the 50° to 75° range, by as much as 
20 dB due to the partial illumination effect of the 50-ns 
pulse. However, Plant's curve at 6 m/s compares very well 


with the (|b la 


rec trans 


2 r 
) curves for the Pierson-Neumann and 


Pierson-Moskowitz models with a maximum deviation of 3 
dB and the (Jb /a 


rec trans 


’) curve for the Elfouhaily Spectrum 


has an overall deviation of 3 dB. 


The results for vertical polarization shown in Fig. 5 are 
also similar to those shown in Fig. 4 in that the simulated 
curves exhibit good agreement with the curve from Long 
within the range of elevation of O°to 40°. 


| == Elfouhaily Spectrum 
Ss L —— Pierson-Moskow itz Spectrum 
-20L | = Pierson-Neumann Spectrum 
=¢- 4FR--Daley 8910 MHz 5 m/s hh (1971) 
=æ- Long 8.8-10 GHz hh (2000) 
E == Long 8910 MHz 5 m/s hh (2000) 
-40- = Plant Show ex X-Band Fixed U = 4 m/s hh (2005) 
| mår Plant Show ex X-Band Fixed U = 6 m/s hh (2005) 


-80+ 


Pye rane! > (4B) 
$ 


-100+ 


0 10 20 30 40 50 60 70 
Elevation Beam Steering Angle (Degrees) 


Fig. 4. Relative Received Power vs. Elevation Beam Steering Angle of a 
horizontally polarized circular dish above a fully developed sea surface with 
a wind speed of 5 m/s 


In examining Fig. 5, it is seen that the larger deviation from 
Long's curve is as much as 7 dB for all of the simulation 
curves, which is due to the response within the enhanced 


b /a 


rec trans 


portion. Again, the ( ’) curves generated from the 


Elfouhaily spectrum appear to be the overall better estimate 
for the higher elevation angles, with a smallest 1-db 
deviation compared to the lower bound and a 5-dB deviation 
compared to the upper bound. 
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O EEDS N a pa 
= Elfouhaily Spectrum 

Ss —— Pierson-Moskow itz Spectrum 
= Pierson-Neumann Spectrum 


-20- | <a Long 8.8-10 GHz vv (2000) 
== Long 8910 MHz vv (2000) 
E =œ- 4FR--Daley 8910 MHz vv (1971) 
-40- |= Plant Show ex X-Band Fixed U = 4 m/s vv (2005) > 


tee Plant Show ex X-Band Fixed U = 6 m/s vv (2005) 


-80- 


2. 

<Irocl rans! > (48) 
È 
o 


-100- 
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Fig. 5. Relative Received Power vs. Elevation Beam Steering Angle of a 
vertically polarized circular dish above a fully developed sea surface with a 
wind speed of 5 m/s 


A special note about the simulated data is that the 
integration error took precedence whenever the grid spacing 
was not small enough (a major cause for the apparent 
oscillations in the curves). Another note is that for smaller 
angles of incidence, a few of the data points are missing, 
which is probably due to the queuing of the parallel 
processor. 


VI. THE GNRCS AS A METRIC FOR SPATIAL COHERENCE 


In Figs. 2 through 5, the relative power received from the 
sea echo, as predicted by the GNRCS, was in substantial 
agreement with experimental baseline values. These curves 
described the relationship between the surface and 
electromagnetic properties within the footprint. However, 
given the size of the footprint and the horizontal correlation 
lengths of sea surface, these curves were expected to 
converge since in most cases the footprint was limited by the 
antenna while the range extent associated with pulsewidth 
was limited by the beam's extent. Hence, the footprint 
consisted of sufficient number of correlation lengths such 


that the GNRCS was essentially identical with o° . 


100000: 


Down Range 4 
Cross Range H 7 
Pulse Length 50 ns 
Pulse Length 500 ns 


ere 


10000; 


Range Extent Compared to E.M. Wavelength, à 


I 
0 10 20 30 40 50 60 70 
Beam Steering Elevation Angle (Degrees) 


Fig. 6. Range Extent vs. Beam Steering Elevation Angles curves for the 
antenna footprint dimensions (cross range and down range) and pulse 
lengths (50 ns and 500 ns) 


Consider Fig. 6 where the cross-range and down-range 
extents, and the range extents associated with two 
pulsewidths are shown (in terms of electromagnetic 
wavelengths) as a function of elevation. It can be seen where 
the footprint is represented by a beam-limited or pulse- 
limited illumination profile. The pulsewidth range extent 
(down-range extent) curves are shown for 500 ns and 50 ns. 
The 500-ns curves show that for nearly all elevation angles 
the footprint is beam limited, representing a fully illuminated 
surface for each pulse. This would mean that for the 500 ns 
pulse the NRCS (o°) is a measure of the coherence effects 
of the entire scattering surface defined by A. However, the 
50-ns curves demonstrate that the footprint is pulse limited. 
Thus, the NRCS values given by Plant would signify that it 
was only partially illuminated, and hence the NRCS could 
only be defined ad hoc by the associated range extent for the 
appropriate angular coverage. Overall, these curves help to 
show that if the number of correlation lengths is sufficient to 
cover the antenna-beam footprint, then o° is all one would 
need to describe the scattering nature of the surface. 
However, if the dimensions of A were small enough to limit 
the number of correlation distances within it, would o° still 
remain a good scattering metric? 


To address this question, consider the product of the gain 
and the back scattering cross section 


G(0,¢)o(S), 


where the o(S) is given by (19). If o(S) is replaced by 


o°, then this product would be essentially linearly 
proportional to the gain. Therefore, the extent to which the 
plots of the two products versus gain differ can be taken as a 
metric of the effects of spatial coherence on the 
backscattering cross section. In other words, if the two 


products are identical, then o(S)=o°. The product 
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G(0,¢)o(S) is computed in Figs. 7 and 8 at a fixed beam- 
center position for a circular dish when the normalized radius 
of the dish (k,a) varied from 50 to 200 for a given 
polarization is plotted. 


The G(0,¢)o(S) curves are shown using o(S) 


obtained from the Elfouhaily Spectrum. For each 
polarization, the G(6,¢)o° curves of Daley and Long 
follow a linear trend. For horizontal polarization in Fig. 7, 


the two baseline curves are separated by 5 dB with a spread 
of 10 to 17 dB on the low end and 23 to 27 dB on the upper 


end. For gain values below 39 dB, the G(@,¢)o(S) curve 
deviates substantially from the linear trend but it is nearly 


linear and well above the baseline curves for gain values 
above 39 dB. 


= Elfouhaily Spectrum (Small Slope 5m/s hh) 
—@ Long 8910 MHz 5 m/s hh (2000) 
© 4FR--Daley 8910 MHz hh 5 m/s (1971) 


36 38 40 42 44 46 
Gain (dB) 


Fig. 7. The product of Gain with a(S ) (o°) vs. Gain for horizontal 


polarization (with a wind speed of 5 m/s) 


For vertical polarization in Fig. 8, the three baseline curves 
are separated by as much as 4 dB with a spread of 13 to 17 
dB on the lower end and 23 to 28 dB on the upper end. The 


G(0,¢)o(S) curve here follows a similar trend as for the 


horizontal polarization. Both of these figures show that as the 
gain increases, or the antenna-footprint dimensions decrease, 


G(0,¢)o° increases linearly. However, the G(0,¢)o(S) 


curve has a non-linear trend that is well above the upper- 
bound baseline curve, which implies that the scattering 
crossing section is gain dependent. Hence, the number of 
horizontal spatial correlation lengths establishes the relative 
strength of the backscattered field. This demonstrates that the 
GNRCS is the best metric for characterizing radar clutter. 


10+----: | = Elfouhaily Spectrum (Small Slope 5m/s vv) 
=e 4FR--Daley 8910 MHz vv 5 m/s (1971) 
—@® Long 8910 MHz 5 m/s vv (2000) 
—@® Long 8.8-10 GHz vv 


36 38 40 42 44 46 
Gain (dB) 


Fig. 8. The product of Gain with a(S) (o°) vs. Gain for vertical 


polarization (with a wind speed of 5 m/s) 


VII. CONCLUSIONS 


The coherence effects on the backscattered field due to a 
random sea surface, as established by the properties of the 
radiating antenna, were investigated using the antenna- 
reciprocity relationship to predict the mono-static signal 
returns from beam-resolved distributed targets for sea-clutter 
analysis. The sea-surface backscatter was examined for a 
fully developed wind-driven sea modeled by several 
empirical wave-number spectra in the presence of a 
uniformly illuminated circular dish at X-band. The relative 
power returns were determined for various antenna 
parameters (height above surface, polarization, and beam 
footprint) for a given sea state. The relative-power curves 
generated by the GNRCS compared very well with responses 
derived from measured data. The deviations obtained from 
the data-derived baseline trends were as small as 3 dB. 
Oddly enough, the simulated responses performed 
unexpectedly well at the high-incidence angles, a major 
indication that rough-surface scattering analysis should 
include the full nature of the incident field as well the surface 
kinematics. Overall, the results show that spatial-coherence 
effects can significantly influence radar-clutter estimates. 


A special note about the simulated data, which was 
previously mentioned, is that the integration error took 
precedence whenever the grid spacing was not small enough, 
which is a major cause for the apparent oscillations in the 
curves. 
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APPENDIX A 


DERIVATION OF THE ANTENNA RESPONSE FROM A FIELD 
SCATTERED BY AN IRREGULAR DIELECTRIC /AIR 
INTERFACE USING A MODIFIED VERSION OF THE SSA 


Fig. A1. Description of the scattering media 


The bounding surface between the dielectric medium and 
free space in Fig. A1 is defined with respect to the reference 


plane H by H —¢(x, y). The free-space and dielectric 
regions correspond, respectively, to z < H -minç(x, y) and 
z2 H-maxç(x,y). A plane waveE,e ™:®-1"7, incident 


from free space together with the outgoing waves, is 
represented by 


= = jk pik, — ky p+ jka (z-2H) 
E(r)=E,e +E, .n€ 


se that Ja’k, E, (kky jee er”), 


(A1) 
In accordance with the Rayleigh hypothesis, the last term, 

represented by the Fourier integral, and referenced to the 

plane z = H , is an approximation to the continuous spectrum 


and is defined as the scattered field, excluding the specular 


-= jk; jk,,(z-2H m š 7 a 
Bea p+ Heal? ) The electric field transmitted into 


term Epe 

the dielectric region is 
(2) = 

E (r) = Eep 


pe IaH fark, pê) (kk, err er, 


e jka H o7 IKu -P-jkzzi(z-H) 


(A2) 
and similarly for the corresponding magnetic fields H(r) 
and H” (r). i, (k,,k,,) and IO (k,,k,,) are electric- 
field scattering amplitudes, respectively, in air and the 


dielectric medium. The following is an outline of the 
procedure used to find an approximate solution to (A1) for 


B(k,,k,) (essentially a Taylor expansion to O|V,c|) and 


relate it to the generalized NRCS, which is one of the 
objectives of the paper under review. 


The symbols indicated here and in the subsequent 
equations are defined as follows: k, is the transverse wave 


number (i.e., the projection of the wave number k on the xy 


plane); k, and k,, are z-components of k , respectively, 
in air and the dielectric medium; k, is the free-space wave 


number; subscript i refers to quantities associated with 
incident fields, subscript t to transverse components; 


and Ky =Z) XK , are unit vectors along the x and 


y axes, p is the position vector in the xy plane, and 


= |A 
Zo = 


& ° 


Matching the six tangential components of the total 
electric and magnetic fields on the media bounding surface 
and, independently, the tangential components of the plane 
wave on the hypothetical planar boundary H , yields the 
following set of integral equations for the scattering 
amplitudes 


E, (k,,k,) -jak K (k k,) K,,(k,,k,) || E, (k;,k,) 
i.) (k ky) ' K,,(k,,k, ) K,, (k,,k,) i.) (kky) 
V, (k k,) 
at 
(A3) 


where the matrix elements K „n (k,,k, ) are 3x3 dyadics that 


depend on wave numbers, the surface slope, and plane-wave 
Fresnel coefficients. Upon applying the Neumann iteration 
procedure to (A3) one finds that to represent the solution to 


within O(|Vc}) the first term suffices. This is equivalent to 
u,(k,,k,) by V,(k,,k,) in (A3), which is 
(it is important to note this approximation does not constrain 
the size of |k,c]) 


approximating 


i, (k, k;) =T? "Ve (k, k,)+ TË "Va (k.k; ); 


(A4) 
where 
=k, € k ke 
T? (k,) = zr K 4 z Ep [aer ok. b 
k,é, +k, “ee k, +k, eoe k,é, +k, ” 
(A5) 
and 
—k,k,, 1 kk. 
T® k z oz ; a a ; 
H ( a) Ere k, +k,, Kafi kê (k,e, +k) ° . 
(A6) 
and 
E E, xz ese) 4 Eo) XZ esy) 
Ve (ha) =7 y epee ft a) a As fl o) 
Qn -(E_. xz, )e 77% 
refr 0 
(A7) 
and 
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p | (20xH Jem astes) +(x Hoy Je! Rus(y) 


-(z, x Hap )e jkanis(x:y) 


(A8) 


Recasting this in spherical coordinates, we denote the 
angles of incidence by 0, (elevation), and ¢, (azimuth) and 
and 8, the electric-field 
scattering amplitude resolved along the spherical unit vectors 
a, and £, assumes the form 


1, (kky) = fd? pe~ 


the scattering angles by a@ 


‘MW (k kap); (A9) 


where 


Y (k,,k;,p) = fa, [Tu E+ ibe E, | + B, eer Eig + Top, Ep } 


(A10) 
and 
Tyo (k, k cos(8-¢,) 
ti? 
a) (ye sma é,-sin’ ate, cosa) 
E cos @, —,/e, -sin Ve, -sina |( eas (O) L o7 Fons " 
eka c(p 


— jk 
+l, cos 6, +4/£, —sin pe ae -e none) 


(A11) 
sin( B-¢ 
Tag (k,, k,,,p)= ( ) e 
(27) erie cos] 
[ 4. =e, =sin? æ cos 9, |(e* eřezsl a 
+e, +6, =sin? a cosg, |T" (0, )(e e~ kus(P) T 
(A12) 
-sin( 8 -¢,) 
T, k,,k,, p = . 
hie (2) (Je, -sma +cosa) 
[1- Je, =sin? æ cos een —efstt)) 
+{1+ E, -sin? «cos, |T" 0, e~ tkuslP) _ pjtous(e) 
eae l i ) (A13) 
-cos(8-¢,) 
T, k,,k;,p = . 
da (22)' (fe, -sin? a + cosa) 
[ Ve, -sin? a -cos (er -e70 ) 
+ é,-sin’@ +cos 6, |r" 0, (e0) _ eas) l 
ee (o ) (A14) 


Next we consider the reception of the scattered field by an 
antenna. From the generalized reciprocity relationship [23] 
we have 


(k,)d?°k,. (A15) 


ee = Qk zzl rad 


b 


: ; ae 2, : 
With this normalization |b,,.| is the received power, and 


the electric far field on transmission is 


tall) ne 


A16 
R (A16) 
with the corresponding radiated power 
p| d?k, 
Fa k : oC Girans A17 
J] i ( ‘) F =] ; oy) 


For compatibility with these definitions E(k, ) represents 
the Fourier Transform of the scattered field at z =0 , i.e., 


eke 


- fE, (p.oje™*a'p 


(A18) 


To use formula (A15) we must express E(k, ) in terms of 


the scattering amplitude E,(k,,k,,)in (A1). For a single 


incident spectral component (i.e., a plane wave) 


S Pai 2 - jk, -p „—-jk,(z-H) 
E; (k i PZ ,) ffa k E i, (k;,k,)e e i (A19) 
For a general illumination we need 
E,(p,0)= (k,,,0,0,)d°k,, 
(27 (A20) 


which is restricted to the plane z=0, the location of the 
antenna. Substituting (A20) in (A18) we get 


3 (k,k, Joi kM 9 Nk)? a? gk dk, 


l, (k, k, Je =" (k, +k, )d’*k,d’k, 


3, (-ki,k, Je ak, 
(A21) 


Substituting the last expression in (A19) we get the 


receiver output in terms of the spatial spectrum: 


=j F, (k 


[Jek 


Cauet e 
2ak,JZ, 


trans “rec 


(A22) 
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A form more suitable for stationary phase evaluation is 
3, (-k,,k,,) in (A22) with (A9) 


obtained by replacing 


_-ifa P| [Ena (k; 


Y (-k),k,,p)e*"e * ®d’k,d’k,. 


trans bec 
2zk,4]Z, 


(A23) 


For k,R>>1 ((A20) is good enough) we can use 
stationary phase, one corresponding to eR the other to 
e *® with stationary points at a=0, B=¢,and 6 =0, 
¢, =@. This translates into k, =k,, = @,k,sin@ and in (A9) 

=-4,; B, = 


—ġ, . As result we get 


-j2kR 
£ R? Fad (k,)- 
2k, JZ, 


(-ak, sin 8, ak, sind, p) 
a 


trans rec 


(A24) 


The matrix elements in (A11), (A12), (A13), and (A14) 
evaluated at the stationary points are: 


2j 


Tya (-k),k,,,2) ~ Gaye 0)sin| k,¢(p)cosé ]; 
To (kikap) =- "(@)sin[ k,o(p) cosa |; 


Tyo, (kiku, P) = Tas, (KK, 2) = 0, 


(A25) 
so that 
Y (—ak, sin 0, ak, sind, p) = 
2J ee ' " 
J ~sin| k,o(p)cos o |{-0,r (A) Eg + G0 (O)E,,\. 
(27) (A26) 


Here as in (A9) E and E, are proportional to F,,,, and 


Fag» the components of Faa (k,) in (A16). We can find the 


proportionality constants by recognizing that (A16) is the 
asymptotic form of the electric field E,,,(R)as k,R—00 


that can be computed from 


1 pe , 
E „a (R) =—— |Ei(k, Je *“"d’k,. (A27) 
J(R)= > Js) 
A stationary phase evaluation gives 
E,.q(R)~ ko paer cos 6E; (0,9) 
rad (22) k,R , 
= JZ Ft (OP) oror, 
° R (A28) 


Writing (A28) in component form 
k? 


2 ai oso O Eo +P Eis), 
WA l ky l ') (A29) 


Fao + DF, radġ 


we obtain the desired proportionality constants: 


_ {Zo 27 F aap E JZ 20F rad 


id ig 7 


jk, cos jk, cos (430) 
Using these in (A26), 
k 0 
¥(-ak, sin, ak, sind, p) VZ sin[ os (P )cos ] 
mk, cosO 
{FOE (8) Foto + $00" (8) Fas}: (A31) 


The final form of the receiver output is obtained by the 
substitution of (A31) into (A24) 


— j2kyR k 0 
GuransOre <= j2nfd? pcos 2h sinf 2A Jos Jp. (99) [- -0,T' (0) Fao + Or" 


coso |fe r (0 )F; 


ee ado 7T" (0) Fy, 


27 eR k cos 
a4 a 0 
jZ jep = 


k, R? 


(A32) 


Qt p, eR k, cosO "Tr(a 0 
AnansD, a ra e ka in[ kos (p) Jose] p “I (2) T al 
2H py _@ MP 
= J] 0a Fe (8-9) Sar (9-9) Fra (0,9) 


APPENDIX B 


BRIEF DERIVATION OF THE ROUGH- SURFACE FIELD 
INTEGRAL EQUATION FOR THE MODIFIED SAA 


A synopsis of the derivation of some of the key formulas 
used to obtain the integral equations given in Appendix A 
and the scattering matrix given in (9) of the paper are 
explained herein. The formalism given in Appendix A will 
be adopted in this document. We wish to emphasize that the 
formulation in our paper is an alternate version to the 
formulation given in [16], although ours renders substantially 
different expressions with unique functional properties. The 
differences in the expressions are primarily due to the 
application of the surface constraints for small slope 


\V.o(x,y)|<«1 and small amplitude |k, (x, y)| <1, and the 


phase reference associated with the mean surface 
displacement (which will be explained in the sequel). The 
methodology presented in [16] and [25], and summarized in 
[26] utilizes the surface constraints simultaneously for the 
SSA. However, the formulation given in Section VII and 
herein, employ these constraints separately. 

Herein, the approximate form of the rough-surface 
scattered field is obtained by a construction that assumes the 
validity of the Rayleigh hypothesis (See Fig. B1). Beginning 
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with the description of the aforementioned construction of 
the electromagnetic fields in both media, we have the 
following. 


(0 


Angle of Incidence 


9) 


Incident Field Scattered Field 


Components Specular 
(Rough Surface) Reflected Field 
Components 
(Mean Level 
Surface) 


Mean Level of Displacement 


Rough-Surface 
profile 


Fig. B1. Description of rough-surface scattering phenomena 
The electric field in the region defined by z < H -ç (x, y) 
(Medium 1), illustrated in Fig. B1, can be written as 
EY (r) = E,e Mee tz 4 Byer TOB) 


pe fJark, B, (k,,k,, ere ae, 


The incident and reflected plane waves, E,e *? "” and 


ei p+ jk, (z-2H) 


E pef , respectively, are characterized at the 


hypothetical plane z = H , and 
Q, (kr) =e" ffa’k, 


scat 


i - jk, -p+ jk, (2-H 
“s (k,,k, Je” AE (B2) 


is the rough-surface scattered field, which is also evaluated at 


the plane z = H . Similarly, the field quantities in Medium 
2 are 


E® (r)=E 
+e" f fd’k, 


e" elke “p— ky (2-H) 


refr 


B® (k,,k, Je Pee), 


(B3) 


= jk, A Hk P— ikon (2-H) 


Here, Eppe” e is a refracted plane wave 


traveling in Medium 2 and the field scattered into Medium 2 
is written as 


BË, (kpe) =e [fak 


pO (kk, ) ePi (H) 


(B4) 


Note that the phase term e/" was included in the 
expressions for the specular fields at z = H to establish the 
phase reference at that location. This permits the reflected 
and refracted plane wave components (as well as 


and E® (k,,k,,)) to be proportionate to the incident field 


via the Fresnel Reflection Coefficients. The purpose of 
eo ikall 


can be explained easily explained in the context of 
transmission-line theory, where the antenna (input) port is 
the reference terminal plane and the load plane located at 
zZ=H. Here, the reflection (transmission) coefficient at 
load plane is referenced to the input reflection coefficient as 


(B1) 


represented in the specular terms of (B1) and (B3). We have 
also assumed that the rough surface scattered fields, as 
expressed using the Fourier expansion theorem, are valid 
within the region (i.e., the Rayleigh Hypothesis) 


H-ming(x, y)<z< H -maxç(x, y). 


The boundary condition for the specularly reflected 
tangential electric-field components at H (i.e, ç(x,y)=0 


) is given as 


(E, + Ee -Eet )% Zo =0. fa) 


The corresponding magnetic-field expressions for both 
media have the given form of 
HO (r) = Hie fie }kaz +H eo “pt jk; (2-2H) 


+e f fd’k H, (k,,k, Je?) 


(B6) 
for Medium 1, and 
BH) (r) = Hy eka g7 IKu P- jka (2-2H) 
- jk; H 2 (2) — jk, -p—jky,(z-H) 
+e" f fd’k HO (k,,k, Je a 


for Medium 2. In the above equations, the plane-wave 
spectral representations for the scattered magnetic fields are 
written as 


(1) aka — jk, p+ jk, (z-H) 
HO, (kpr) =e" [[d°k HI, (k,,k, Je?" (aes 


and 


HO), (kr) = gan [Ja?k a? (kk, je" enn, (B9) 


Similarly, the boundary condition for the tangential 
components for the magnetic plane wave fields at the mean 
level z= H is also given as 


(H, +H.) -Hp )x% = 0. (B10) 


In examining the phase term in (B2), (B4), (B8), and (B9), 
e/" is treated as an arbitrary constant submitted in the 
expressions to give rise to a certain symmetry. It can also be 
said that each member of the spectral ensemble for the 
scattered fields is phase referenced to the planez=H. 
Again, this step yields a symmetric set of solutions for the 


scattering amplitudes E®®” (k,,k,) andH{”(k,,k,), 
which vanish when¢(x,y)=0. Note that the scale of 


|s(x,y)| has yet to be specified; only the Rayleigh 


Hypothesis was enforced. The scattering amplitudes have a 
set of impedance relationships, which were derived from 
Maxwell's Equations and are given as 


(k, —k,Z, )x 


k,,k 


t? >ti 


7 
a ( 


) = ou H, (k,,k;) 


(B11) 
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and 


(k, + k, 2, )x 


i) (K,,k,,) =ou, HY (k,,K, ). 


(B12) 


Having defined the necessary parameters and field 


quantities in each medium, the solutions to the scattering 
amplitude forE,(k,,k,), and hence H,(k,,k,) can be 
determined by using (B1), (B3), (B6), and (B7) at 


z=H-¢(x,y) with the boundary conditions for the 


tangential field components, which are given by (B5), (B10), 
(B11), and (B12). These steps are shown as follows. The 
tangential electric-field components must be continuous at 


z =H-¢(x,y). Therefore, the boundary condition imposed 
on the surface is represented by the following expression, 


[(E; X Vo jern +E, XV G les? 


+{Ja’k, 


= (Exes XV, je RE + fJa’k, 


ry ~ JK, p- jk,s(x.y 

i(k k; )x ye Ve ee) 

q(2 — Fk, pt jkao(%Y 
a?) (k ky Jx ye errr), 


(B13) 


for the electric field and 


[(v xH, yeast) + (v, x Hen jer? | eP 


+[fd’k,v, xH, (k,, k, ) e0 


2 (v, x He jg terie) + ffa "kV, x HO (k k; jeer een 
(B14) 
for the magnetic field. The outward normal is defined as 
+V 
v, =- A G (B15) 
0 2 
1+ (V.S) 


Next, the cancellation of ,/1+V,.¢-V,¢ on both sides of 
(B13) and (B14) is performed and the normal v, is now set 
to-v, =z,+V,¢. For convenience, the integration variable 
was changed from k, tok,. A further modification is 
required to exploit a particular feature of atx y) in the 
expressions above, in that 


“jer 


= [ak vox 


a ' ~ — jk, 
d, (k, k, )x ve EA 


a, (kky Je? + ffd’k;z, xE, (k,,k, Je? 


(B16) 


and 


[fer 


= [Jark,V,sx 


i (kik, Jx pe? 


a.) (kk, Je? + [fd7k\z, x EO (k, k, )e™? 


(B17) 


are added and subtracted from both sides of (B13), 
respectively. Equation (B13) is then re-expressed as 


E, xV, ese) Ea XV gasy) | po ikue 

(E: x7.) (En Yo) 
+f fa’k; 
-Í fd’k;V,gx 


-f| fa’k;z, x 


=(E 


B, (kik, xpe (eM) 1) 


a (kik, Je? 


i (ksk, Jem” 


ee (B18) 
xV, jew pt iky2iS (XY) 


refr 


+f fark, 
—[[a°k Vox 
-f fd’k;z, x 


i) (kky xpe? (en -1) 


P (kpk, Je? 


JO (i,k, Je, 


Correspondingly, the magnetic field has a similar form to 
the above 


[(%, xH, Jett) A (¥, Hen jer) | P 
+ Sjak [kx (kk, )]e e (ese = 
070 


1 
k,Zo 
1 


[fa°k.v,cx[ kB, (kk, )Je™* 


[Jark.2, x 


= (v xH, Jepretan 
0 refr 


= [k'xE,(k;,ky ) Je” 


0 


(B19) 


; [Jark.y, x k; x BY) (kku )] ee (eres) -1) 


0 


1 


+ 


0 


[larkv.cx Ík; xE® (kok; ) ]e** 


070 


1 


ffd’kiz, x| k, xB? (kky) |e. 


k,Z, 


Ho 
E 


where k, is the free space wave number and Z, = . The 


vector wave numbers k` and k, are defined as 
k =k, -k,z,; 


; , 7 B20 
k, =k, -k,,Zo- i i 


The impedance relationships given in (B11) and (B12) 
were used to obtain (B19). Note that the terms containing the 
jkaag (xY 


+ jk,¢(x,y) 


factors e —1 and e )_1 can be shown to be 


negligibly small compared 10 O(|V,5 (x, yJ). We now 


ik,- 
L grer 


(27) 


multiply both sides of (B18) with and integrate 


over the entire xy plane. Since, 


jo d’p=(2r)' 6(k, -k,) (B21) 
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we obtain 


~ jkag(x,y) Fi ~ this (X.Y) 
Ja pe? (E; xv, jem + (Even Vo )e aa 


(Ey x Py Jetse y) 
+) fa°k, fa? pe" i }e i, (k,,k,,)x¥, Ga y) 1) 
t +V cox, (k,,k,,) 
— fark: fa? pels)? O (kku) x, (ers 1-1) 
t +V ox © (kk) 
=(22) {29 x 3 (2) (k,,k,;)—Z, x 3, (koku )} 


(B22) 
for (B18). A similar step performed on (B19) yields 


(¥, xH. „Je es x,y) + (Vv, % Ha )e -jkag(xy) 


jere 


; )] ( eis) -1) 

(ki ky) 

: )|( eikiss(Y) -1) 
(ks) 

i (k,,k,) |}. 


(B23) 


With these given expressions, the final solution for the 
vector scattering amplitudes can be obtained by using a 
series of tedious algebraic steps (which will not be shown 
here). The final result is a solution to the Fredholm Integral 
equation of the second kind 


k,Z,(d 2 pelt ele 


+f fa?k; fa? pe 


-f fa?k; fa? pe“ K)o 


= (27) fz, x[kx 


i, (k,,k,) fer | K,,(k,k,) K,,(k,k,) |] B,(k,,k,) 
EP (k,,k,) i K,,(k,,k, ) K.,, (k,,k, ) a” (k,,k,, ) 
A V (k ,k;) 
V, (k,,k;)[ 
(B24) 


where the first-order solutions for the scattering amplitudes 
are 


ti 


ti 


w (B25) 


ka) | _ 
ka) | 
V; ( P-V, (kok,) 
O vy (koki) 


V, (k,, a)+T 


and driving functions V, ,, 


k Eyer 
k 
(k,, 


k, ) are defined as 


(E, xý, Jese) + (Even xv, )e -jkas(x,y) 


(Ey x V ) etas) 


k,-k, 
et) e 


Ve (k,,k,) z 


| 


(B26) 


1 
(27)’ 
and 


(¥, xH, Je etast) +(¥, Hen Je -jkas(x,y) 


-7 xH, Je jkonis(xy) 


k,Z, 


_ _o“o 


(B27) 


Recall that under the  small-slope constraint, 
\V.s(xy)/<<1,%, ~1; and as a result, (B26) and (B27) 
reduce to (A7) and (A8) (in Appendix A), respectively. In 
examining (B22) and (B23), one should be reminded that 


there is no constraint on the scale of [kos (x, y) 


, and thus, 


the continuity of the tangential electromagnetic field 
components hold at z = H —¢(x, y) as long as the Rayleigh 


Hypothesis is valid. Another consideration is that ¢(x,y) is 


completely independent of H , and thus the scattered field 
lacks any phase continuity related to the difference between 


¢(x,y) andH . Therefore, any level shift perturbation of 
¢(x,y) (say by some constant +h,) would be completely 


absorbed into H , and the formulation given in Appendix A, 
and herein would still apply. 


Jimmy O. Alatishe was born in Washington, DC, in 1975. He received his 
B.S., M.S., and Ph.D degress in Electrical Engineering in 1998, 2000, and 
2012, respectively, from The George Washington University in Washington, 
DC. His primary interest are in the areas of rough surface scattering and 
microwave circuits. He has been employed as a research engineer by the 
Radar Division at The Naval Research Laboratory in Washington, DC, 
Since 2000. His current efforts include research on the development of 
Radar recevier and transmitter modules, and advanced signal processing 
techniques. 


Wasyl Wasylkiwskyj (F'97) received his BEE degree from the City 
University of New York in 1957, and the MS and Ph.D. degrees in electrical 
engineering from the Polytechnic University in 1965 and 1968, respectively. 
In 1984, after a 15 year career in private industry, where his last position 
was vice president and general manager of Physical Dynamics, Inc. in La 
Jolla California, he joined the George Washington University, Washington, 
DC, as Professor of Engineering and Applied Science. His research 
experience in the academic environment as well as a consultant to industry 
covers a broad spectrum of applied electromagnetics. This includes 
microwave components and techniques, phased-Array Antennas, EM 
propagation and scattering, radar cross-section modeling as well as 
modeling geophysical and oceanographic electromagnetic phenomena. In 
1997 he was elected Fellow of IEEE for “original contributions to 
electromagnetics, propagation and scattering and to fundamental 
understanding of phased-Array Antennas”. His relativelty recent research 
was in direction finding algorithms with particular emphasis on the effects 
of mutual coupling on algorithm performance. The focus of his current 
research interests is in antennas for geophysical applications and scattering 
from random surfaces. 


16 


